Enhanced carrier scattering rates in dilute magnetic semiconductors with correlated impurities 
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In III-V dilute magnetic semiconductors (DMSs) such as Gai-ajMn^ As, the impurity positions tend to be 
correlated, which can drastically affect the electronic transport properties of these materials. Within the memory 
function formalism we have derived a general expression for the current relaxation kernel in spin and charge 
disordered media and have calculated spin and charge scattering rates in the weak-disorder limit. Using a 
simple model for magnetic impurity clustering, we find a significant enhancement of the charge scattering. The 
enhancement is sensitive to cluster parameters and may be controllable through post-growth annealing. 



PACS numbers: 72.80Ey, 78.30Ly 
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The perspective of utilizing charge and spin of the electrons 
for new electronic device applications has generated tremen- 
dous interest in the field of spintronics 1 1]. A unique combi- 
nation of magnetic and semiconducting properties makes di- 
lute magnetic semiconductors (DMSs) very attractive for var- 
ious spintronics applications |2]. Among the family of DMSs, 
much attention has been paid to Gai-^Mn^As since the dis- 
covery of its relatively high ferromagnetic transition tempera- 
ture 01, with a current record of T c — 159 K 

In Ga 1 _ :c Mn ;E As, unlike in II- VI DMSs, the magnetic ions 
in substitutional positions act as acceptors delivering one hole 
per ion. All Gai-^Mn^As samples are, however, heavily 
compensated, with hole concentrations much less than x. 
This signals the presence of substantial amounts of donor de- 
fects like arsenic antisites AsQa or interstitial manganese ions 
Miii generated during low temperature molecular beam epi- 
taxial growth 0]. The magnetic and transport properties of 
Ga 1 _ a; Mn a ;As depend not only on the manganese fraction x 
but are extremely sensitive to detailed growth conditions |5], 
as well as to temperature and speed of post-growth annealing 
eILZIIMI- This sensitivity points to the crucial role played by 
the defects and their configuration, and has stimulated intense 
research on the structure of defects and their influence on the 
magnetic and transport properties of DMSs |9, 10]. 

Most theoretical models for transport in DMSs assume ran- 
dom defect distributions. However, the presence of both pos- 
itively and negatively charged defects results in a correlation 
of their positions. Indeed, Timm et al. |9] found in the limit 
of thermal equilibrium that, driven by Coulomb attraction, the 
defects tend to form clusters. The main effect of clustering is 
ionic screening of the disorder Coulomb potential, which has 
been shown to be necessary to correctly reproduce the band 
gap, metal-insulator transition and shape of the magnetization 
curve 1 9]. 

In this paper, we demonstrate that the correlation of defect 
positions can have a dramatic effect on electronic transport in 
DMSs: the conductivity of Gai^Mn-r As is strongly modified 
through a momentum dependent impurity structure factor. We 
will show that the clustering significantly increases the charge 
scattering relaxation rate. At the same time, positional cor- 
relation taken alone is not sufficient to affect spin scattering: 
orientational correlation of spin scatterers is also necessary. 

To describe the transport properties of DMSs we have em- 
ployed the memory function formalism fill fl2l 11311 . The 



central point of this approach is the calculation of the cur- 
rent relaxation kernel (or memory function), whose imaginary 
part can be associated with the Drude relaxation rate. To get 
an expression for the memory function in spin- and charge- 
disordered media we have used the equation of motion ap- 
proach 1141 M 511 (details of the derivation will be presented 
elsewhere). Here, we are particularly interested in a param- 
agnetic system in the long-wavelength limit, the case relevant 
for studying the conductivity in DMSs above T c . In this case 
the memory function is obtained as 
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where n is the carrier concentration and the operators p* 1 are 
defined through a four-component (/z = 1, +, — , z) charge 
and spin density vector, 
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Here, c^ 1 is defined via the Pauli matrices, where a 1 is the 
2x2 unit matrix, a ± = (a x ±ia v )/2, and Xp^p" (k, w) are the 
associated charge- and spin-density response functions. The 
superscript c in Eq. Q refers to a clean (defect-free) system. 

The presence of impurities, including their correlations, en- 
ters in Eq. Q through the expression k)£/„(k))jy m . 
The angular brackets indicate a thermodynamical average 
with respect to a magnetic subsystem Hamiltonian H m . We 
assume H m to be a sum of individual spin contributions cor- 
responding to uncorrelated and noninteracting localized spins. 
The disorder scattering potential is described by the four- 
component impurity charge- and spin-density operator 
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where the summation is performed over all defect positions. 
In order to separate the effect of the impurity structure fac- 
tor from other effects of clustering like ionic screening we 
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consider only one type of defects, Mn 2+ ions in cation sub- 
stitutional positions. All defects in Eq. produce thus the 
same charge potential and carry localized spins. The latter 
are treated as quantum mechanical operators coupled to the 
band carriers via a contact Heisenberg interaction resulting in 
a momentum-independent exchange constant J. In our calcu- 
lations we neglect the time evolution of localized spin opera- 
tors, thus assuming that the carriers move through an ensem- 
ble of frozen spins. For the charge component E/i(k) we take 
a Coulomb potential screened with the host material dielec- 
tric constant. Effects of ionic screening are disregarded, and 
screening by the electron liquid is absorbed in the response 
functions. 

Let us now consider the correlated product of two compo- 
nents of the disorder potential in Eq. Q. First, we separate 
the same-ion (J = j') and pair (j ^ f) contributions: 
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where m is the impurity concentration, and in the second term 
the average of the product becomes a product of averages for 
noninteracting spins. Random impurities are taken into ac- 
count through the first term on the right-hand side of Eq. ©, 
while the second term contributes only if there are correla- 
tions in the impurity positions. For spin scattering, however, 
spatial correlations are not sufficient for the pair term to sur- 
vive. Indeed, if any of the indices \i, v corresponds to a spin 
component (e.g., [i — z), then the second term in Eq. is 
proportional to the average spin and vanishes if (S z ) = 0, re- 
gardless of spatial correlations. The presence of two sources 
of randomness is a characteristic feature of spin scattering. In 
terms of scattering, localized spins are correlated if they ex- 
hibit both positional and orientational correlations. Note that 
the presence of a macroscopic magnetization is not necessary 
for spins to be correlated. What counts in the scattering is 
the short-range orientational correlation that might be present 
even in a macroscopically paramagnetic system. 
For the charge-scattering term (/i, v — 1), we have 
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with the structure factor 



a 



Vx ^ 



(5) 



(6) 



Mj'' 



where x is the molar fraction of magnetic ions in the sample 
and £l is the elementary cell volume. Let us introduce a pair 
distribution function P(R), normalized as 



1 J P(R)dR = x, 



(7) 



which describes the probability to find another magnetic ion 
at a distance R from a given ion. We approximate Eq. as 

S(k) w 1 + — / P(R) cos (k • R) dR. (8) 



For a random impurity distribution one has P(R) = x, and 
the second term in Eq. (jSJl vanishes. The structure factor 5(k) 
is then equal to 1, which implies a contribution only of the 
same-ion term in Eq. @. 

To study the effect of correlations in the defect positions 
(higher probability to find magnetic impurities close to each 
other), we employ a simple model expression for a pair distri- 
bution function -P(R), assuming it to be a piecewise constant, 
spherically symmetrical function of the form 



P(R) = 



R < R c , 

R c < R < Rd 

R > Rd ■ 
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The first region corresponds to a cluster of radius R c with ef- 
fective impurity concentration x c > x, the second region is 
a depletion layer with x c < x, necessary to preserve the av- 
erage impurity concentration in the sample. The width of the 
depletion layer is determined by the normalization condition 
(0. The impurity concentration within the depletion layer has 
a minor impact on the final results, and to keep things simple 
we fix it to Xd = x/2. 

The remaining two parameters, R c and x c , describe the 
cluster structure and are in general independent. We can re- 
late them, however, if we fix the average number N of the 
impurity ions within the cluster: 



'-x P . = N. 
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This seems reasonable for modelling the effect of annealing 
on low-temperature grown DMS samples, where one may as- 
sume that the total number of Mn 2+ substitutional ions within 
the cluster is conserved while the cluster size (and density) 
may vary. In the following we will use N = 10, consistent 
with the results of Monte-Carlo calculations |9]. 

Our model pair correlation function (|9} yields the following 
momentum-dependent impurity structure factor: 
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Fig-fflshows S(k) for two different values of the cluster radius 
R c , which in principle can be controlled by annealing. As ex- 
pected, the structure factor oscillates with decreasing ampli- 
tude, with a larger first maximum for smaller cluster size. 

Eq. Q contains the set of full system charge- and spin- 
density response functions and, strictly speaking, should be 
calculated by iterations. This approach was realized in 1 16] to 
study a spin-independent system close to the metal insulator 
threshold. In our case, however, we assume that the disorder is 
weak enough so we can approximate Eq. Q by expanding to 
second order in the disorder potential U (k), and thus replace 
Xpfp" (k, uj) by its clean system counterpart Xpv- P " (k, oj). 

An accurate description of carrier-mediated ferromagne- 
tism H7II and optical response IU8I1 of Gai-^Mn^As would 
require taking the true multiband structure of the material into 
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FIG. 1 : Momentum-dependent impurity structure factor i 1 H for dif- 
ferent values of cluster radius R c . Number of ions within the cluster 
N — 10, average impurity concentration in the sample x = 0.05. 
Inset: schematic plot of the pair distribution function {5). 



FIG. 2: Frequency dependence of charge i 141 and spin 1151 scatter- 
ing relaxation rates, for impurity clusters with R c — 10A, x c = 0.1. 
Inset: illustration of the frequency dependence of the cluster en- 
hancement, see text. 



account. In this paper, however, we concentrate primarily on 
effects arising from the disorder configuration, putting less 
emphasis on the details of the band structure. In the evalua- 
tion of the response functions Xp^p" (k, oj) we thus work with 
a simple parabolic band. 

Furthermore, we use the random phase approximation 
(RPA) to account for electron-electron interactions. For sim- 
plicity we limit ourselves to a static RPA and express the 
density-density response function of the interacting system as 



Xnn(q,w) = 



xo(q,^) 
£RPA(q, o) ' 



(12) 



where xo(qj w ) 1S the non-interacting density-density re- 
sponse function of the clean system, i.e., the Lindhard func- 
tion, and £RPA(q, 0) is the static RPA dielectric function 1 19]. 
The spin response functions in the paramagnetic state are not 
affected by electron-electron interactions on the RPA level, 
and they can also be expressed in terms of the Lindhard func- 
tion xo(q,w). 

Under the above approximations Eq. Q can be directly 
evaluated and the memory function becomes the sum of spin 
and charge contributions 
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Here Sun — 5/2 is the localized spin of magnetic impurities, 
and the common prefactor is given by A = (m / n)V 2 /&ir 2 m. 
The imaginary parts of Eqs. dl4> and ill 5b represent the energy 
dependent charge- and spin-scattering contributions to the 
Drude relaxation rate. Generally speaking, within our model 
the relaxation times T n i s \ (oj, q) are also momentum depen- 
dent. In the present paper, however, we consider only a long- 
wavelength limit, setting q — > 0. Note that the momentum- 
dependent impurity structure factor S(k) does not appear in 
the spin term (1151 . 

We have evaluated the scattering rates, Eqs. (I14> and d!5i . 
for the case of Gao.95Mno.05As with a hole concentration of 
p = 0.5 hole per magnetic ion. Other parameters used are: 
heavy hole effective mass rn = 0.5 m , dielectric constant 
e = 13, and exchange constant VJ = 55 meVnm 3 , which 
corresponds to the widely used DMS p-d exchange constant 
NnB= 1.2eVfii. 

In Fig. 13 we plot the frequency dependence of the imagi- 
nary part of Eqs. J14I and J15I . Both spin and charge relax- 
ation rates demonstrate the frequency dependencies that one 
might expect for momentum-independent scattering and for 
Coulomb scattering. In agreement with earlier estimations by 
golden rule 1 20], the charge scattering in dc limit dominates 
the spin scattering. We would like to point out, however, that 
in our calculations (as well as in Ref. 1 20]) we did not take 
into account the effect of ionic screening shown ||9j] to reduce 
significantly the charge disorder potential. In any case, even 
in the present model the spin scattering contribution is clearly 
not negligible, and reaches the same order of magnitude as 
charge scattering for higher frequencies. Both contributions 
should therefore be taken into account simultaneously. 

Correlation in impurity positions results in a significant in- 
crease of the charge scattering contribution at low frequency, 
while for higher u this enhancement decreases. The origin of 
this effect is illustrated in the inset in Fig.|2] The presence of 
impurity clusters selects excitations within a finite momentum 
window defined by the impurity structure factor. On the other 
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FIG. 3: Relaxation rate enhancement <16l due to correlation in im- 
purity positions as a function of cluster size. The average number of 
magnetic ions within the cluster is fixed to iV = 10. 



hand, the region of one-particle excitations is given by the 
imaginary part of the Lindhard function \o (shaded region in 
Fig.|2ji. The maximum enhancement of the relaxation rate cor- 
responds to maximum overlap of the selection window with 
the one-particle spectrum, which happens at low frequency. 
For higher ui the window falls outside of the available excita- 
tion spectrum, reducing the relaxation rate enhancement. 
In Fig.|5]we plot the cluster enhancement factor 



is fixed to N — 10, and the cluster radius R c is related to 
concentration of magnetic ions within the cluster x c through 
Eq. J10I . The enhancement is strongest (up to a factor of 3) for 
low frequency, and is quite sensitive to cluster configuration. 
The latter will be sensitive to post-growth annealing, which 
is widely used to increase T c in Gai-^Mn^As samples. The 
possible modification of transport properties described here 
should therefore be taken into account along with other effects 
of annealing like decrease of the number of interstitial Mn 
ions and increase of carrier concentration. 

To summarize, using the memory function formalism we 
have considered transport in charge and spin disordered media 
with potential application for DMSs, with particular emphasis 
on non-randomness of impurity positions in Gai-^Mn-cAs. 
We have shown that positional correlations alone of the mag- 
netic impurities do not affect spin scattering: orientational 
correlations are also necessary. For charge scattering, im- 
purity clustering gives rise to a momentum-dependent impu- 
rity structure factor which substantially modifies the transport 
properties of the material, typically leading to 100% enhance- 
ments of relaxation rates. These results should give valu- 
able insights into the effects of annealing on low temperature 
grown DMS samples. 

Finally, the discussion in this paper was limited to DMS in 
the paramagnetic state. However, it is well known that mag- 
netic ordering in DMS can have a dramatic influence on trans- 
port properties |21]. Our approach should be well suited to 
study these effects. 



(t*)-i + (t*) 



R\-l 
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as a function of cluster configuration. The enhancement factor 
is defined as the ratio of the total (charge plus spin) relaxation 
rates for correlated and random impurity distributions. Recall 
that the average number of magnetic ions within the cluster 



Acknowledgments 



This work was supported by DOE Grant No. DE-FG02- 
05ER46213. 



[1] G. A. Prinz, Science 282, 1660 (1998). 
[2] H. Ohno, Science 281, 951 (1998). 

[3] K. Edmonds, P. Boguslawski, K. Wang, R. Campion, N. Farley, 

B. Gallagher, C. Foxon, M. Sawicki, T. Dietl, M. Nardelli, and 

J. Bernholc, Phys. Rev. Lett. 92, 037201 (2004). 
[4] H. Munekata, H. Ohno, S. von Molnar, A. Segmuller, L. L. 

Chang, and L. Esaki, Phys. Rev. Lett. 63, 1849 (1989). 
[5] H. Shimizu, T. Hayashi, T. Nishinaga, and M. Tanaka, Appl. 

Phys. Lett. 74, 398 (1999). 
[6] T. Hayashi, Y. Hashimoto, S. Katsumoto, and Y. Iye, Appl. 

Phys. Lett. 78, 1691 (2001). 
[7] S. J. Potashnik, K. C. Ku, S. H. Chun, J. J. Berry, N. Samarth, 

and P. Schiffer, App. Phys. Lett 79, 1495 (2001). 
[8] K. M. Yu, W. Walukiewicz, T. Wojtowicz, I. Kuryliszyn, X. 

Liu, Y. Sasaki, and J. K. Furdyna, Phys. Rev. B 65, 201303(R) 

(2002). 

[9] C. Timm, F. Schafer and F. von Oppen, Phys. Rev. Lett. 89, 
137201 (2002). 

[10] X. Y. Cui, J. E. Medvedeva, B. Delley, A. J. Freeman, N. New- 
man, and C. Stampfl, Phys. Rev. Lett. 95, 256404 (2005). 



[11] W. Gotze, Phil. Mag. B 43, 219 (1981). 

[12] D. Belitz and S. Das Sarma, Phys. Rev. B 34, 8264 (1986). 

[13] C. A. Ullrich and G. Vignale, Phys. Rev. B 65, 245102 (2002); 

Phys. Rev. B 70, 239903(E) (2004). 
[14] W. Gotze, P. Wolfle, Phys. Rev. B 6, 1226 (1972). 
[15] G. F. Giuliani and G. Vignale, Quantum Theory of the Electron 

Liquid (Cambridge University Press, Cambridge, 2005). 
[16] A. Gold and W. Gotze, Phys. Rev. B 33, 2495 (1986). 
[17] T. Dietl, H. Ohno and F. Matsukura, Phys. Rev. B 63, 195205 

(2001). 

[18] J. Sinova, T. Jungwirth, S.-R. Eric Yang, J. Kucera, and A. H. 

MacDonald, Phys. Rev. B 66, 041202(R) (2002). 
[19] G. D. Mahan, Many-Particle Physics, 3rd ed. (Plenum, New 

York, 2000). 

[20] T. Jungwirth, M. Abolfath, J. Sinova, J. Kucera, and A.H. Mac- 
Donald, Appl. Phys. Lett. 81 4029 (2002). 

[21] T. Omiya, F. Matsukara, T. Dietl, Y. Ohno, T. Sakon, M. Mo- 
tokawa, and H. Ohno, Physica E 7, 976 (2000). 



